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1.  Introduction 


Recently,  combined  methods  have  been  introduced  for  the  determination 
of  the  deflections  of  the  vertical.  Those  methods  may,  first  of  all,  incorporate 
the  following  data:  satellite  derived  potential  coefficients,  mean  gravity  anomalies 
and  point  gravity  observations.  The  idea  is  that  the  inner  zone  is  determined  in 
detail  from  the  point  anomalies,  for  instance,  by  using  the  method  of  least-squares 
collocation,  while  the  effect  of  an  intermediate  zone  is  calculated  with  the  Vening 
Meinesz'  formula  using  mean  gravity  anomalies,  and  finally,  the  remote  zone 
contribution  is  represented  by  a spherical  hannonic  expansion. 

In  this  paper  we  are  going  to  study  the  errors  of  the  combined  method. 
The  investigation  is  mainly  restricted  to  the  errors  of  the  intermediate  and  remote 
zones.  The  error  contribution  of  the  inner  zone  differs  being  much  dependent  on 
the  method  of  computation  and  density  and  quality  of  the  data.  For  this  zone,  we 
will  adopt  some  representative  values  found  in  the  literature. 


2.  Computational  Method 

We  briefly  summarize  the  combined  method  presented  in  Lachapelle 
(1977).  The  components  4 and  7)  are  considered  to  consist  of  three  subcomponents 
(^0,  ii,  nofVifVz)  such  that: 


(1) 


4 - 4o  + 4i  + 

T?  = T?o+  7?1  + 7)2 


and 


Each  of  the  subcomponents  is  the  contribution  of  the  deflections  from  a specific 
zone  around  the  point  of  computation  (see  Figure  1).  Exceptions  are  4o  and  tjo  , 
which  also  serve  as  "reference  field"  in  di  and  da  . 


Figure  1.  Subdivision  of  the  Surface  Around  the  Point  of  Computation 
into  Inner  Zone  (d^).  Intermediate  Zone  (da)  and  Remote 
Zone  (do ). 
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The  values  of  ^ o and  77 o are  computed  from  the  fully  normalized  potential  coeff- 
icients Jna  and  Kna  in  the  following  Way  ( Lachapelle,  ibid.,  p.  3): 


inJ  ■i(r)  V [j"»  S- 

^ n =2  m = 0 


where  na.x  is  the  maximum  degree  of  expansion,  and; 


j-  Rna  1 ^ r 

D(p 

'^Sn.(<P,X)J  ~ 1 

r Rn«(*P  » ^ r 

{-  ^ } = 1 

- sin  mX . 
CO&  mX‘ 


„a  (sincp) 


m 


cos(p 


ra/r  = ratio  between  the  radius  (ra)  of  the  internal  sphere 
(to  which  Jna  and  Kn#  are  referred)  and  the  radius 
(r)  of  the  point  of  computation. 


Remark  1.  Formula  (2)  is  a slight  generalization  01  Lachapelle's  formula  for  points 
at  an  arbitrary  height  above  the  sphere  to  which  the  coefficients  J„t  and  Kna  refer. 
For  low  degree  expansions  at  the  surface  of  the  earth  rs/r  = 1 is  a good  approxi- 
mation. 

The  values  of  4a  and  rja  are  obtained  by  applying  Vening  Meinesz' 
integral  formula  to  gravity  anomalies  that  are  formed  by  subtracting  from  mean 
terrestrial  free-air  anomalies  (Ag)  the  contribution  (Ags)  implied  by  the  potential 
coefficients  used  in  computing  4o  and  T)o  . This  computation,  which  is  carried 
out  for  the  zone  Os  with  il>i^  ^ t can  be  expressed  in  the  following  way: 


where  S{il))  is  Stokes'  function,  do  is  an  elemental  area,  and  a is  the  azimuth 
from  the  point  of  computation  to  a current  point. 
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The  classical  Vening  Meinesz'  formula  is  valid  strictly  only  for  a 
sphere.  However,  terrestrial  data  Ag  may  be  used  provided  that  either  they 
are  continued  analytically  to  the  internal  sphere  or  that  a terrain  correction 
(Molodenskii  term)  is  applied  to  Ag  (Heiskanen  and  Moritz,  1967,  pp.  315  and 
313,  respectively). 

Lachapelle  (ibid.)  used  l°x  1°  mean  gravity  anomalies  for  these 
computations,  which  were  extended  to  = 8°.  ibi  varied  between  0? 7 and 
1°5,  depending  on  the  density  of  the  point  gravity  data  in  the  inner  zone. 

The  values  of  ^ i and  T)i  may  be  obtained  in  different  ways; 
Lachapelle  (ibid.)  used  least  squares  collocation.  Such  an  approach  has  the 
advantage  that  there  is  no  difficulty  to  compute  the  effect  of  the  innermost  zone 
and  we  can  easily  incorporate  heterogeneous  data  in  the  computations.  As 
pointed  out  by  Tscheming  (1974)  a considerable  gain  in  accuracy  will  be  achieved 
if  we  include  astrogeodetic  deflection  of  the  vertical  in  the  set  of  observations. 

On  the  other  hand,  we  have  to  limit  the  number  of  observations  to  a few  hundred 
in  order  to  keep  the  computer  time  at  a reasonable  level.  Another  difficulty  is 
the  instability  that  occurs  if  two  observations  are  located  close  to  each  other. 

The  original  Vening  Meinesz'  formula  does  not  suffer  from  these  limitations. 


3.  Error  Analysis 


We  now  study  the  error  propagation  due  to  the  data  and  other  sources. 
It  is  assumed  that  the  numerical  integration  of  Vening  Meinesz'  formula  is  per- 
formed with  such  an  accuracy  that  the  integration  errors  can  be  neglected. 
Furthermore,  the  model  errors  caused  by  disregarding  the  flattening  of  the 
earth  (in  the  order  of  3 x 10”^)  are  not  considert*d. 

The  gravity  data  are  assumed  to  be  corrected  for  the  terrain  effect  in 
Vening  Meinesz'  formula.  If  this  correction  is  omitted,  an  additional  error  in  the 
order  of  0'.'2,  is  introduced  (Moritz,  1966;  Dimitrijevich,  1972). 

The  effect  of  the  earth’s  atmosphere  can  be  estimated  in  the  following 
way  (Moritz,  1974).  A constant  distribution  of  air  is  assumed  above  the  reference 
ellipsoid.  This  homogeneous  mass  has  no  effect  on  the  deflection.  Next,  the 
fictitious  density  of  the  atmosphere  in  ide  the  topography  (alxive  the  referent'e 
ellipsoid)  is  subtracted  from  the  density  of  the  tojMigraphy.  (This  modified  density 
of  the  topography  can  be  used  for  a simultaneous  correction  for  terrain  and 
atmosphere.)  We  obtain  the  following  contribution  from  each  compartment  to  the 
correction  for  the  atmosphere* 


atm.  corr.  = - terrain  corr.  x 


density  of  atmoS[)here 
density  of  to))ogni))hy 
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As  the  terrain  correction  is  in  the  order  0'.'2  it  is  obvious  that  the  atmosphere  has 
no  practical  effect  on  the  deflections  of  the  vertical. 


Finally,  we  regard  the  earth  constants  as  known  (without  errors). 


3.1  Errors  from  the  Remote  Zone 

The  error  contribution  of  the  remote  zone  is  dependent  on  the  spherical 
distance  (is  (see  Figure  1).  A formula  for  the  RMS  influence  of  his  zone  on  the 
total  deflection  of  the  vertical: 


(4)  e = V 

is  given  in  Heiskanen  and  Moritz  (1967,  p.  262).  However,  as  pointed  out  by 
deWitte(  1966)  and  Hagiwara  (1972),  Molodenskii's  truncation  coefficients,  Q^,  in 
this  formula  need  to  be  modified.  A generalization  is  obtained  by  substituting  c„ 

by  (ra/r)^^"*^^  c„  (cf.  Remark  1).  With  these  modifications,  the  mean  square 
contribution  of  the  remote  zone  becomes- 


(5) 

o» 

CD 

to 

II 

n = 3 

where 

c„ 

anomaly  degree  variances  at  the  internal  (Bjerhammar) 

sphere 

Qn(^2)  = 

Qn  (0a)  + S(02)  P„  (COS  02)  sin  02  / n(n+  1) 

S((i)s)  = 

Stokes'  function 

il>a 

geocentric  distance  of  truncation 

As  the  spherical  harmonic  expansion  for  a given  set  of  coefficients  can  be  regarded 
as  a reference  field  in  the  inner  and  intermediate  zones,  errors  in  the  coefficients 
will  influence  the  remote  zone  contributions  only.  Thus,  the  errors  of  the  remote 
zones  are  of  two  kinds:  the  potential  coefficient  errors  and  the  truncation  errors. 
From  (2)  and  (5)  we  obtain  the  following  mean  square  propagation  of  the  potential 
coefficient  errors: 


^ ffla 


(6) 


66 


a(n+3) 


0,1 


= Y n (n+ 1)  Q ® ((ta)  (^)  dc„ 


where  dCn  is  the  mean  square  error  of  Cn.  This  error  can  be  determined  in  the 
following  way. 

Let  us  expand  the  gravity  anomaly  (Ag*)  at  the  internal  sphere  (re  = r) 
into  a series  of  spherical  harmonics  to  degree  nj,*,  [cf.  (2)]: 

II  max 

(7)  = Y Ag*n 

n = 2 

where 

n 

Ag*  = G(n-  1)  Y [Jn*.  RnmCP,  >•  ) + K„.  S„„  (<P  , X )] 

D=  0 


Thus,  the  coefficient  errors  6 and  SKna  are  propagated  according  to: 

n 


6Agt  = G(n-l)y  R„.(<P,X)+  6Kn«S„,  ((p.X)] 

m=  0 


and  the  global  mean  square  error  of  Cn  is  finally  given  by: 


(8)  dc„=  M{6Agt®]  = G^(n- 1)^  Y (6jf.  + 

m — 0 


where  M {x]  is  the  global  average  of  x: 

= s JJ " 

a 


In  (8)  we  have  used  the  orthogonality  property  of  the  spherical  harmonics  over  a 
sphere.  For  a different  estimate  of  6cn,  see  Rapp  (1973).  In  Table  1 we  give 
the  mean  square  errors  6cn  computed  from  the  errors  of  the  GEM  7 coefficients 
(Wagner,  1976). 
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Table  1.  Mean  Square  Errors  of  the  Degree  Variances 
for  GEM  7 from  Wagner  (1976,  Table  27). 
Fully  Normalized  Harmonics.  G = 980  gal. 


n 

6a®  10^®* 

6c„  [mgal^Jt 

2 

35 

0.0000 

3 

357 

0.0014 

4 

219 

0.0019 

5 

938 

0.0144 

6 

715 

0.0172 

7 

2103 

0.0727 

8 

1710 

0.0805 

9 

3368 

0.2070 

10 

3191 

0.2503 

11 

5990 

0.5753 

12 

4924 

0.5722 

13 

7865 

1.0877 

14 

7247 

1.  1762 

15 

11271 

2.1216 

16 

11357 

2.4541 

* 6cr^  = y + 6k®  ) 

-J 

D =:  0 


t 6c„  = G^  (n-  1)^ 


Finally,  by  inserting  these  values  for  6cn  into  formula  (6),  the  potential  coefficient 
error  for  n„ax  = 16  was  determined.  The  ratio  re/r  was  set  equal  to  1 in  this  low 
degree  expansion.  The  result  is  depicted  in  Figure  2. 


The  other  error  source  of  the  remote  zone,  the  error  due  to  the 
■ n>-  . . n of  the  spherical  harmonic  expansion,  is  also  given  by  formula  (5). 
mean  square  value  of  the  truncation  error  is: 


The 


(9) 


60 


0.2 


4G" 


n(n+  1)  qI 


(^3)  (— ) C, 


, + l 
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For  the  numerical  computation  of  this  error,  we  use  the  degree  variances  c 
of  Tscheming  and  Rapp  (1974): 


(10)  c„  = A (n- l)/(n- 2)(n+ 24)  , n^S 

s = {re/rf  = 0.999617 
A = 425.28  mgal^ 


The  convergence  of  (9)  is  very  slow  for  <1)2  close  to  0 while  for  larger  angles, 
say  5°,  it  converges  well.  Formula  (9)  is  illustrated  in  Figure  2.  The 
excellent  subroutine  of  Paul  (1973)  was  used  for  a rapid  determination  of  , 
and  the  series  was  truncated  at  n = 2000.  A remarkable  minimum  (0'.'06)  of  the 
RMS  error  is  obtained  for  02  = 40°. 


3.2  Errors  Due  to  Lack  of  More  Detailed  Data  in  Oa 


The  application  of  Vening  Meinesz'  formula  (3)  for  the  determination 
of  la  and  rja  requires,  theoretically,  that  Ag  is  known  at  each  point  in  (Ja.  In 
practice  we  limit  ourselves  to  using  mean  gravity  anomalies  in  this  area.  We 
shall  now  estimate  the  error  due  to  the  lack  of  more  detailed  gravity  material  in 
CTa  • The  data  is  assumed  to  be  located  on  a mean  earth  sphere  of  radius  r. 

Let  us  expand  Ag  into  Laplace  harmonics: 

00 

Ag  = y Ag„ 

r.  = 0 


The  mean  gravity  field  (^)  is  related  to  Ag  according  to: 


^ (y) 


CO 

JJ  B (yx)  Ag  do  = y Agn 


n=2 


where  B(yx)  is  the  averaging  operator  and  /3n  are  its  eigen  values.  If  we 
approximate  each  block  of  mean  anomalies  with  a circular  cap  with  equal  area  we 
obtain  (see  Meissl,  1971  , p.  24): 
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Figure  2.  The  RMS  errors  of  the  Total  Deflection  of  the  Vertical  Due  to 
Potential  Coefficient  Errors  and  Truncation  Error. 


= [Pn-iJCOSV'o)-  Pn  + i (COS  O )] /(2  D + 1)(  1- COS  o) 


where  i/Jq  is  the  spherical  radius  of  the  cap.  If  v is  the  block  size,  then  is 
given  by : 


il)o  = i>  sin  u/-n 


Now  the  error  at  each  point  when  representing  Ag  by  is: 


6g=Ag-Ag  = y (l-i9„)Ag„ 


Following  the  derivation  of  formula  (5)  in  Heiskanen  and  Moritz  (1967,  pp.  261- 
262),  we  obtain  the  following  error  propagation  of  6g  outside  the  spherical 
distance  i/). 

The  error  of  each  component  ^ V becomes; 


COSO 


and  the  total  RMS  error  is  given  by: 


66  = M { 64^  + 6n^]  = 


CO  CO 

= A ^ Q.(i/')Qn>(<^)  m{ 

40  -i  ->  t 


d 6g,  ^6gn'  1 ^ „ 36g  1 _ 


30  ^O  cos  O SX  9X 


A y QU<i>)  n(n+l)  M {6g^  ] 
4G 
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where 


6g„  = (1-  ^„)  Ag„ 

and 

M{6g!]=  (l-/3„f  M {Ag^  } =(l-i3„f 


Hence,  by  using  Vening  Meinesz'  formula  outside  the  spherical  distance  4) , the 
following  mean  square  error  (due  to  neglecting  more  detaUed  data)  is  committed: 


= «())  n(n+l)(l-i3„)2 

4G  -• 

n = 2 


Cn 


Finally,  the  corresponding  contribution  from  the  zone  Q2  is  given  by: 

00 

(12)  = [Q^d^i)- Qn(</J2)]  n(n+ l)(l-i3^) 

4G  -■ 

n = 3 


As  approaches  0 for  large  n this  formula  suffers  from  the  same  slow  conver- 
gence as  (9)  for  small  angles  Ji)i . For  4°  (12)  was  expanded  to  n = 5000. 

For  larger  , n = 2000  was  found  to  be  a sufficient  degree  of  truncation  (1^=  1°). 
The  results  are  shown  in  Figure  3. 


3.3  Errors  Due  to  Inaccurate  Gravity  Material  in 


A constant  error  in  Ag  will  not  affect  the  deflections  of  the  vertical, 
because  there  is  no  zero-order  term  present  in  ^ and  t?.  Hence,  the  only  error 
source  that  has  to  be  considered  is  the  random  errors  due  to  insufficient  gravity 
data  within  each  block  to  assure  an  accurate  determination  of  the  mean  value. 
The  propagation  of  this  error  is  given  by  (3).  We  obtain  approximately: 


(13) 


66 


a, a - 


(4ttG)‘ 


I [JJ 


d^ 


d a 


CTi 
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Figure  3.  The  RMS  Error  of  theTotal  Deflection  of  the  Vertical  Due  to 
Lack  of  More  Detailed  Gravity  Material  than  1°  x 1° 
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where  m i is  the  standard  error  of  ^ i and  a i is  the  area  of  block  i.  The 
evaluation  of  (13)  requires  that  m i is  known  for  all  blocks  used  in  the  deter- 
mination of  is  and  77s.  It  will  therefore  be  very  much  dependent  on  the  quality 
of  the  mean  anomalies  and  will  vary  for  each  point  of  computation.  In  Table  2, 
we  report  the  computations  in  four  specific  points.  The  computations  are  based 
on  the  l°x  1°  anomaly  information  described  in  Rapp  (1977).  For  unknown 
blocks  m 1 is  set  to  30  mgal.  In  this  table  we  also  show  the  error  estimates  of 
Groten  and  Moritz  (1964).  By  assuming  uniform  errors  of  all  ^ blocks,  they 
arrived  at  the  following  formula: 


(14) 


^®a,a  ~ s J(<i'2)] 

8^(Gr) 


where 


J(4») 


I 


Ip 


sin  ^ d 


S = 0.027  r^  (l°x  1°  blocks,  1 gravity  profile  inside  each  block). 
Formula  (14)  seems  to  give  a reasonable  approximation  of  662,2  for  Ipi  ' 2^ , 
specially  in  areas  with  a poor  gravity  material  (points  2 and  4), 

In  Table  2 we  also  report  the  sum  of  squares  of  the  errors  considered 
above  (the  error  of  the  inner  zone  is  not  included).  From  these  results  we  conclude 
that  the  RMS  errors  are  significantly  decreasing  with  an  increasing  angle  ^s  all  the 
way  to  30°. 

3.4  An  Extended  View 


On  the  basis  of  the  previous  computation  results,  it  is  also  of  interest  to 
study  the  errors  when  ips  is  extended  beyond  30°.  However,  such  an  extension 
would  include  large  areas  where  no  l°x  1°  gravity  material  exists  today.  In  these 
computations,  we  will  therefore  assume  that  we  have  such  a material  with  a uniform 
error  of  all  blocks.  Thus,  we  can  apply  formula  (14)  for  the  computation  of  663,2. 
The  resulting  RMS  errors  are  shown  in  Figure  4. 

From  Figures  2 and  4 we  conclude  that  a dominating  error  source  is  the 
truncation  of  the  spherical  harmonic  expansion.  A significant  minimum  of  the  RMS 
error  is  obtained  for  ~ 40°.  It  is  unreasonable  to  extend  liis  beyond  this  minimum. 


-13- 


3.5  The  Inner  Zone  Error  and  the  Total  Error 


The  error  of  the  inner  zone  depends  on  the  quality  and  distribution  of 
the  observations  in  this  area.  To  some  extent  it  will  also  vary  with  the  method 
used  in  the  computations.  By  usii^  least  squares  collocation  with  approximately 
200  observations,  Lachapelle  (ibid.)  estimated  the  errors  of  4i  and  r)i  to  about 
1'.'3.  These  errors  were  rather  independent  of  the  cap  size  , which  varied 
between  0?7  and  1°5  in  the  computations.  Hence,  by  using  this  technique,  the 
total  RMS  error  of  0 will  be  in  the  order  of  2:'2  for  4>i  - 1°  and  ibs  = 40°. 

Kearsley  (1976)  applied  the  original  Vening  Meinesz'  formula  with 
the  well-known  subdivision  into  Rice  Rings.  The  inner  zone  errors  of  and  T)i 
were  about  01'3  for  ^ 1° . Thus,  a total  RMS  error  of  0 in  the  order  of  1'.'3 
(for  (l)a  = 40°)  seems  possible  to  achieve  in  an  intensive  determination  of  the 
inner  zone  contribution.  Even  though  the  same  accuracy  might  be  possible  to 
obtain  by  using  least  squares  collocation,  there  are  practical  limits  of  the  number 
of  observations  in  this  method  (Lachapelle,  ibid. , p.  6).  Here  we  refer  to 
6^=  677=  0'.'3  as  the  accuracy  of  Vening  Meinesz'  formula.  Some  total  RMS 
errors  based  on  the  above  estimates  of  6 are  given  in  Table  3. 


Table  3.  Total  RMS  Errors  of  the  Deflections  of  the  Vertical 
for  = 1° 


m 

0 

0 

iH 

1 

0 

0 

1 

30*" 

0 

0 

Collocation  69^=  ll'84 

2';88 

2'.'57 

2'.’34 

2'.' 24 

2';22 

Vening  Meinesz  60^  = OL'42 

2 '.'2  6 

i:'85 

1'.'50 

1'.'35 

1'.'32 

Lachapelle  (ibid.)  compared  predicted  deflections  of  the  vertical  at 
some  169  astrogeodetic  stations  with  the  astrogeodetic  deflection  components. 
The  RMS  difference  between  the  total  deflections  was  2:'12  . This  result  agrees 
fairly  well  with  Table  3. 


4.  Summary  and  Conclusions 

In  this  paper  we  have  studied  possible  error  sources  for  the 
determination  of  the  deflections  of  the  vertical  by  using  a combined  method.  The 
investigation  deals  mainly  with  the  errors  generated  in  the  intermediate  and 
remote  zone,  while  for  the  inner  zone,  some  error  figures  have  been  adopted 
from  earlier  studies.  The  numerical  results  are  shown  in  Tables  2-3  and 
Figures  2-4. 
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We  assume  that  the  only  contribution  to  4 and  r}  from  the  remote  zone 
is  a spherical  harmonic  expansion  to  degree  16.  The  RMS  truncation  error  of 
such  a series  was  found  considerably  sensitive  to  changes  in  the  distance  of 
truncation,  <^2  (see  Figure  2).  A local  minimum  (660,2  = O'.'OG)  is  obtained  for 
i^)2  = 40'’. 

In  the  intermediate  zone  the  gravity  field  was  assumed  to  be  repre- 
sented by  l°x  1°  mean  anomalies.  The  error  due  to  lack  of  a more  detailed 
gravity  material  is  depicted  in  Figure  3.  This  error  source  can  be  diminished 
by  using  a more  detailed  subdivision  of  the  blocks  between  the  distances  1°  and 
10°  from  the  computation  pxjint. 

The  error  prqjagation  of  the  mean  anomaly  errors  (662^2)  is  shown 
in  Table  2 for  four  selected  points.  A reasonable  approximation  of  these  errors 
is  obtained  by  using  formula  (14)  according  to  Moritz  and  Groten. 

i The  RMS  sum  of  the  errors  of  the  intermediate  and  remote  zone  is 

given  in  Figure  4.  A significant  gain  in  accuracy  is  achieved  by  extending  ((12  to 
at  least  30°.  A minimum  is  obtained  for  >lis  = 40°.  Finally,  by  adding  a repre- 
sentative error  for  the  inner  zone  computation,  the  total  RMS  error  was  estimated 
(Table  3).  If  the  inner  zone  is  determined  by  the  method  of  least  squares  colloca- 
tion, the  total  RMS  error  will  hardly  be  less  than  21'2  = 1°).  Substituting  this 

method  with  an  accurate  version  of  Vening  Meinesz'  formula,  the  final  error 
might  decrease  to  11'3.  Further  improvements  are  expected  for  a refined  sub- 
division of  the  intermediate  zone  blocks.  In  no  case  should  the  truncation  distance 
(ips)  exceed  40°. 
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